Properties of the Radio-Emitting Gas Around SgrA* 
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We show that the radial profiles of the temperature and density of the electrons as well as the 
magnetic field strength around the massive black hole at the Galactic center, Sgr A* , may be 
constrained directly from existing radio data without any need to make prior assumptions about 
the dynamics of the emitting gas. The observed spectrum and wavelength-dependent angular size 
of Sgr A* indicate that the synchrotron emission originates from an optically-thick plasma of quasi- 
thermal electrons. We find that the electron temperature rises above the virial temperature within 
tens of Schwarzschild radii from the black hole, suggesting that the emitting plasma may be outflow- 
ing. Constraints on the electron density profile are derived from polarization measurements. Our 
best-fit results differ from expectations based on existing theoretical models. However, these models 
cannot be ruled out as of yet due to uncertainties in the source size measurements. Our constraints 
could tighten considerably with future improvements in the size determination and simultaneous 
polarization measurements at multiple wavelengths. 



I. INTRODUCTION 

The supermassive black hole at the Galactic center, 
Sgr A* , occupies the largest angle on the sky among all 
known black holes. Its extended image provides an excel- 
lent opportunity to study the physics of low-luminosity 
accretion flows. 

The bolometric luminosity of Sgr A* ~ 10'^^ erg s~^ is 
~ 8.5 orders of magnitude smaller than the Eddington 
limit for its black hole mass of ~ 4 x 10^ Mq. Over the 
past decade various theoretical models have been pro- 
posed to explain the low luminosity of Sgr A* despite 
the large gas reservoir from stellar winds in its vicin- 
ity. Among the early models invoked was an Advection 
Dominated Accretion Flow (ADAF) involving hot pro- 
tons and cold electrons with a low radiative efficiency at 
the Bondi accretion rate of ~ \Q~^Mq yr~^ [21]. Sub- 
sequently, the detection of linear polarization was used 
to set an upper limit on the electron density near the 
black hole, which ruled out the original ADAF proposal 
Pi [l2 1 22 1 and favored shallower density profiles with a 
lower accretion rate such as in a Convection Dominated 
Accretion Flow (CDAF) Later variants of the ADAF 
model allowed for outflows, namely mass loss from the in- 
flowing gas Most recently, an improved Radiatively 
Inefficient Accretion Flow (RIAF) model was proposed, 
involving substantial mass loss (although the outflowing 
mass is ignored in calculating the radio emission), a non- 
thermal component of electrons, and different electron 
and proton temperatures. Other models associated the 
radio emission with a jet [301] or a compact torus near the 
black hole [H. 

In parallel to these modelling developments, the data 
on Sgr A* has improved dramatically over the past few 
years. The latest observations include new determina- 
tions of the size, spectral luminosity, polarization and 
rotation measure of the source as a function of wave- 



length [U, OJ, 127|, 129|. With the rich data set that is 
now available, it is timely to remove any theoretical prej- 
udice and ask: what does the data alone tell us about 
the properties of the radiating gas? In addressing this 
minimal question here, we deviate from past practice of 
modelers who made assumptions about the dynamics of 
the accreting gas before interpreting the observational 
data on Sgr A* . We avoid dynamical assumptions and 
attempt to constrain the properties of the radio-emitting 
gas directly from the data itself. 

As discussed in detail in § II below, the measurements 
of the source size at different radio wavelengths provide 
crucial constraints on the properties of the gas surround- 
ing Sgr A* . Current size measurements are unfortu- 
nately subject to large error bars, which in turn imply 
large uncertainties in the inferred gas properties. Our 
analysis provides an estimate of the spatial dependence 
of gas properties based on current measurements, adopt- 
ing a frequency dependent size r (x v°' with a = 1 ± 0.3 
l4 \2i\ . Our methodology demonstrates how more ac- 
curate measurements may be used to obtain better con- 
straints with no model-dependent assumptions about the 
dynamics of the gas. 

In the different subsections of § II we apply our ap- 
proach to various aspects of the data on Sgr A* that are 
currently available. We compare our results to previous 
work in § III. Finally, § IV summarizes our main conclu- 
sions. 



II. EMPIRICAL CONSTRAINTS 
II. 1. Radio Spectrum and Size: Data 

The spectral luminosity of SgrA* is time dependent. 
As illustrated in Fig. [TJ the observed specific luminosity 
per unit frequency v at the brighter emission epochs 
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[3, 34 1 is well described by a power-law form in the fre- 
quency range of 3-1000 GHz, 



vL^ = 1.7 X lO^^'i^i^^ erg s'^ 



(1) 



where vn = {v/lQ^^ Hz). Since the flux measured at the 
times used for size determination is close to the brighter 
emission values, we will use this power-law index in our 
phenomenological discussion. 

We note that at low frequencies < 3 GHz the observed 
flux somewhat exceeds the flux given by Eq. ([1]) [equiva- 
lently, fitting a power-law only to the data at < 10 GHz 
would result in a fit that under-predicts the flux at higher 
frequencies (as indicated by the dashed line in Fig. [T]), 
i.e. in a "sub- mm excess"]. This deviation is of no sig- 
nificance to the analysis below, which focuses on higher 
frequencies, > 10 GHz. It does imply, however, that in 
applying our simple analytic results to the lowest ob- 
served frequencies, some minor quantitative modifica- 
tions would need to be introduced. As explained below, 
we argue that the radio flux is dominated at different 
frequencies by plasma located at different radii. The de- 
viation from Eq. ([T]) below 3 GHz implies therefore that 
the gas temperature and magnetic field strength at large 
radii, > lO^**'^ cm, differ slightly from those obtained us- 
ing our simple power-law scalings, which are based on 
Eq. HD. 
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FIG. 1: Specific luminosity of SgrA* from simultaneous multi- 
frequency data of Falcke et al. [3, filled circles], and Zhao et 
al. [SJi up/down triangles representing the flux at different 
times, close to the times of maximum/minimum in the 1 mm 
flux]. Diamonds denote flux measurements by Zylka et al. 
[H, 666 GHz] and Serabyn et al. [23, 850 GHz]. Squares 
denote the flux densities derived from the data used for size 
determination, given in Krichbaum et al. [H, [3 (for 1.4 mm 
and 3.4 mm) and in Shen et al. [l^] (for 7 mm). The solid line 
shows the power-law relation vL,, = 1.7x 10'''*i^ii*erg s~^, and 
the dashed line is 0.7x 10^*!^ii^erg s~^. The assumed distance 
to Sgr A* is 8kpc. 



As demonstrated in § II. 2 below, measurements of the 
size of Sgr A* at various radio wavelengt hs p rovide im- 
portant constraints on the emitting gas 0, [ij, 123] . For a 
black hole mass of M = 4 X IG^Mq 0, the radial scale 
is set by the Schwarzschild radius of Rg = 1.2 x 10"^^ cm, 
which corresponds to an angle of 0.01 mas on the sky 
at our adopted distance of 8 kpc. Table I presents the 
latest data from Krichbaum et al. [14| and compares the 
inferred brightness temperature T;, to the virial temper- 
ature [iH]- The brightness temperature 37) at a radius 
r from Sgr A* where the emissivity at an observed fre- 
quency V peaks, is defined through the relation 



vL^ ^ Anr^iiy) x fg x -^Tb{iy) 



(2) 



where the geometric coefficient /g < 1 is the ratio be- 
tween the emission surface area and the area of a sphere 
of radius r {fg — 1 for a sphere and fg = 0.5 for a 
two-sided disk with the same radius). We conservatively 
assume that the surface area scales as since this pro- 
vides the lowest brightness temperature (which, as we 
will show, is already above the virial temperature at large 
radii). We define the virial temperature Ty by equating 
the thermal kinetic energy of the plasma to half of the 
gravitational potential energy per proton. 



3 GMm^ 



(3) 



where rup is the proton mass. (Note that the "escape 
temperature" at which the thermal kinetic energy of the 
plasma exceeds its gravitational potential energy is 2Ty.) 
Using Eq. (P) and M = 4 x 10'^ Mq we obtain 

n = 7.6(2/5)-V-i-6 



ri3 MeV, 



(4) 



T, = g.Sr"' MeV, 



(5) 



where = (r/10 cm). Table I shows that Tf, must 
be close to at = 1, and that there is preliminary 
evidence that Ti, /Ty increases with radius. 

Shen et al. 27| infer a power-law dependence of the 
intrinsic size of Sgr A* on wavelength A of r cx A'^, with 
P — 1.09 ± 0.33. For our phenomenological analysis we 
use 



-1 n —i/a 
ri3 = l.Oi^^;^ , l/ii 

with a = 13^^ « 1 ± 0.3, implying 
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(6) 



(7) 



II. 2. Radio Spectrum and Size: Implications 

The observed radio luminosity originates most li keh 
from synchrotron emission by relativistic electrons [12| , 
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TABLE L Measured size and brightness temperature of 
Sgr A* . The listed radii are related to the FWHM intrinsic 
sizes of Krichbaum et al. [IJ] by r = 0.5 x FWHM/0.76 (see 
footnote [36]). 
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1.4 
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2.1 
0.88 
0.43 
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The frequency dependence of source size implies that the 
emission cannot originate from an opticaUy-thin plasma. 
To see this, consider the electrons at r dominating the 
emission at a frequency i^{r). If the optical depth is 
small, then these electrons would produce a spectrum 
vLi, oc v'^^^ at < i^{r), close to the observed spectrum 
at these frequencies. This implies that the emission at 
v < v^r) would be dominated by electrons at r, which 
is inconsistent with the frequency dependence of source 
size. We therefore conclude that the optical depth for 
synchrotron self-absorption satisfies T^[iy{r),r] > 1. For 
a similar reason, the characteristic synchrotron emission 
frequency Vcif) = (7e )('3-S/27''mec) of the electrons at r 
[dominating the emission at a frequency v{r)] must sat- 
isfy lydf) ~ v{r). If VciT) ^ v{r) then for Tij[v{r),r\ ^ 
1 these electrons would produce a flux vL^ oc at 
v > inconsistent with the observed spectrum, and 

for Tij[i'(r)^r\ = 1 these electrons would produce a flux 
vLn cx v'^^^ at v > dominating the emission at 

v > i/(r) in conflict with the frequency dependence of 
source size. 

We therefore conclude that radiation at different radii 
r is dominated by electrons with I'df) ~ fir) and that 
T^[i'{r),r] > 1. Next, we argue that the electron energy 
distribution may be characterized by a single energy or 
an effective temperature Tg. That is, we show that the 
energy distribution of electrons cannot be highly non- 
thermal. Consider, for example, a power-law distribu- 
tion of electron energies, dne/dje oc J^^. Such a distri- 
bution would be consistent with observations provided 
that T„[i'{r),r] « 1 (rather than T^[i^{r),r] 3> 1), since 
otherwise the flux emitted by electrons at r would ex- 
tend beyond iy{r) as vL^, oc v"^^^, exceeding the observed 
flux. For Ty[j/(r),r] « 1, vL^, oc i/(3-p)/2 at > i'{r), 
and the value of p is constrained by the ratio between 
the far- infrared luminositj^ ~ 3 x 10'^^ erg/s at a fre- 
quency ~ 10^^ Hz [13, llj, and the radio luminosity, 
~ 5 X 10^^ erg/s at ~ 10^^ Hz, to be p > 4.3. This large 
power-law index implies that only a small fraction of the 
total energy can be carried by electrons of energy exceed- 
ing that of the electrons dominating the radio emission. 
Moreover, we will show in S II. 3 that the extension of 



such a power-law to electron energies much below that of 
the electrons dominating the radio emission would imply 
a very large rotation measure, inconsistent with observa- 
tions |3a|. These constraints are satisfied by recent RIAF 
models^ol [3^ which associate only a small fraction of 
the total electron energy with a power-law component. 
However, our simple analysis shows that if thermal emis- 
sion at different frequencies originates at different radii 
to account for the observed spectrum and size measure- 
ments, then there is no need for an ad-hoc non-thermal 
component. 

Since the emission originates from an optically thick 
plasma, the characteristic temperature (energy) Tg of the 
electrons dominating the radiation is 



1^1.6a-2 
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The electron temperature has to satisfy 



MeV. 



(8) 



Mr) = 12 



Tejr) 



X 0.3 



eB{r) 
2'innpC 



y{r), (9) 



where nif, is the electron mass and we used the relation 
(7g) — 12\Te{r) / nieC^Y which angular brackets denote 
an average over a relativistic Maxwellian of temperature 
Tg. Based on Eqs. ([6]) and ^ this requirement imphes 



Te{r) 



Tv{r) 



B 



(10) 



In deriving this result we have not used the inferred 
source size but only the fact that it is frequency depen- 
dent. Using the size estimate in Eq. ^ and Tg « Tf,, we 
then get 



B 



-4:2a 



13 



G. 



(11) 



The z/(r) oc r " scaling of the observed radia- 
tion frequency on emission radius implies ~ — 
MeV. In order for the synchrotron model to 
hold down to 1 GHz, the electrons must remain rela- 
tivistic, i.e. the condition 5 x lO-^^-'^-^+'^Z") > 1 must 
hold, implying l/a<lorQ;>l. It therefore ap- 
pears that the value 1/a ~ 1 of Shen et al. 27] is pre- 
ferred over the alternative suggestion for higher values 
1/a « 1.5 ±0.2 

The results in Eqs. (O, ([8]) and (fTTj) have several im- 
portant implications. First, Te is close to Ty at ris ~ 1, 
and Te/Ty increases with radius approximately as r^/^ 
for /? ~ 1.1 [3^. This implies that the gas cannot be 
confined to a thin disk and the flow geometry must be 
quasi-spherical. Moreover the radio-emitting gas is not 
likely to be flowing in but rather flowing out since its 
thermal kinetic energy exceeds the gravitational binding 
energy beyond a radius of a few tens of Schwarzschild 
radii (ri3 > 2). Our conclusions would only be strength- 
ened if the emitting plasma follows a jet geometry for 



4 



which the surface area scales as with S < 2 [see the 
discussion following Eq. ([2])]. 

Finally, we note that the plasma under consideration 
is collisionless as the Coulomb collision time is much 
longer than the dynamical time of the gas 27r yVYGM — 

0.9 X W^rH^s. However, collective plasma effects should 
operate [26| since the inverse of the plasma frequency 
or electron gyro-frequency are much shorter than the dy- 
namical time. We use the term "temperature" in our dis- 
cussion to characterize the typical electron energy even 
if the electron distribution function happens to be non- 
Maxwellian. 



II. 3. Density constraints: Opacity, Rotation 
Measure and Circular Polarization 

The optical depth to synchrotron self-absorption is 
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5/0 f ^-l^/.^.-3-0.4a g ^^-3 



,B = 2.0 X 10^(2/g)"V(r)r 
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and 



5.6 X lQ^{2fg)-^T{r)r 



13 



(12) 



(13) 



(14) 



For awl the optical depth increases with radius, r oc 
TLer^ . In order to ensure that rir) > 1, it is sufficient to 
require that this condition will hold at ri3 = 1, implying 
ne(ri3 = 1) > lO^cni^^. 

The relativistic rotation measure of a fluid of electrons 
with a thermal Lorentz factor 7e and density threaded 
by a coherent magnetic field B, is given by RRM = 8 x 
105(ne/cm-3)(B/G)(r/pc)(7-2)rad m-2 As long 

as B is coherent, we may express this rotation measure 
in terms of r as 



RRM 



2.5 X 103(2/g)r(r)r 



,2-3. 6a 
13 

9-7. 4a , 



•ad 



5x 10^(2/g)V6r?3"' "rad 



(15) 
(16) 



where we substituted (7-2) « [771^0^ /Tf. For a k. 1.1, 
the rotation measure scales as RRM oc and may be 
either decreasing or increasing with r. In the latter case, 
the rotation measure is dominated at all frequencies by 
the same outermost electron shell, while in the former 
case it is dominated by electrons near the radius where 
radiation is emitted (and so it is expected to be larger 
for higher frequencies or smaller radii). The observed 
rotation measure of ~ 6 x 10^ rad/m^ at a frequency of 
~ 2 X 10^^ Hz [Ref. [3 and references therein], implies 
for a coherent B-field that 



< 10^(2/,, 



(17) 



Marrone et al. [l3| report measurements at 2.3 x 10^^ Hz 
and 3.5 x 10^^ Hz. While there is an indication that 
the rotation measure is higher at the higher frequency 
by a factor of few (see their Table I) the observations at 
the two frequencies are not simultaneous, and since the 
source is variable the differences may be due to variabil- 
ity. If the rotation measure differences are real and not 
due to the temporal variability, Eq. ()16|) requires 



ne(ri3 = 1) ~ 10^(2/3) 



"^cm ^, 



(18) 



with the density decreasing with r at least as steeply as 
^-9+7. 4a ^ For this density (and fg ^ 1/2), the mag- 
netic field and thermal energy densities are comparable 
at ri3 = 1 and T(ri3 = 1) ~ 10^. A turbulent magnetic 
field would generate a random walk in the net rotation 
measure and so the electron density inferred from the 
RRM observations would increase by the square-root of 
the number of field reversals (coherent B patches) along 
the region where the RRM originates. The large lin- 
ear polarization observed at frequencies v > lO^GHz for 
Sgr A* implies that the inferred magnetic field is not 
highly tangled in the innermost region. The low level 
of linear polarization at lower frequencies is consistent 
with the notion that the emission at different frequencies 
originates from different radii. 

For a « 1 the electrons become mildly relativistic at 
r ~ 10^^ cm, where the emission is predicted to peak 
around ~ 1 GHz. This may account for the circular 
polarization observed at these low frequencies [3, [3, [i^l . 

Finally, we note that (as mentioned in § II. 2) the 
observed RRM excludes a power-law extension of the 
electron energy distribution, due/d'ye oc 7"^ with p > 
4.3, down to energies significantly lower than Tg. The 
contribution of lower energy electrons to the rotation 
measure is proportional to rif,/^^ oc 'ye{dne/dje)/je o^ 
7~P~^. Using Eq. and denoting by 7™ the min- 



imum electron Lorentz factor we have RRM « 2.5 x 
10^T{Te/jm'mec'^)~^~'^ rad/m^, which implies for p > 4.3 
that the rotation measure would exceed the observed 
value of ~ 5 X lO^rad/m^ for 'jmmec'^ /Tf. < 1/2. 

II. 4. Equipartition and Entropy 

The equipartition ratio between the magnetic energy 
density and the thermal energy density of the electrons 
scales as 



and the entropy scales as 



^5 oc . 



T! r 
— oc 

Up n 



4.8a-6 



(19) 



(20) 



Requiring uniform entropy and equipartition fraction 
gives a = 16/14.8 = 1.08 which is surprisingly within 
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the range inferred by Shen et al. This special 

value yields the scalings Tg oc r^°'^^, Ue oc r"*^ ®^, and 
B cx r"" ''^*. Substituting these power-law scalings in Eq. 

implies that the rotation measure RRM is nearly in- 
dependent of emission radius or observed frequency. This 
result can be tested by future observations that would 
monitor the time dependence of RRM at different fre- 
quencies [3| • A similar time dependence at different fre- 
quencies would imply that the rotation measure is dom- 
inated by a common outer shell. 



III. COMPARISON WITH EARLIER WORK 

The RIAF models generically predict that the radio 
emission is dominated by thermal electrons near the in- 
nermost stable circular orbit (ISCO) in the central region 
of the accretion disk (sij. The radio spectrum produced 
in this model by the thermal electrons is inconsistent with 
the observed spectrum, which is well described by the 
power-law form in Eq. ([T]). A power-law electron compo- 
nent is therefore added to the thermal RIAF component 
in an ad-hoc manner, where the spectral index and nor- 
malization of the power-law component are tailored to 
match the low frequency radio data [2^, Is^]. We have 
pointed out in § II. 2 that the size measurements indi- 
cate that the source size is frequency dependent, implying 
that the radiation is dominated at different frequencies by 
electrons at different locations. This in turn suggests that 
the observed spectrum reflects the spatial dependence of 
electron temperature rather than the energy distribution 
of the electrons at a single radius. 

Since the emission of radiation in RIAF models is dom- 
inated by the innermost region of the disk, at r ~ 3i?s — 
3.6 X 10^^ cm, the observed size is dominated in these 
models by foreground interstellar scattering at all fre- 
quencies. The addition of a power-law electron compo- 
nent to these models increases slightly the intrinsic source 
size [33| . but does not change the requirement that the 
measured source size will be dominated at all frequen- 
cies by interstellar medium scattering. These results ap- 
pear to be at odds with the latest size measurements 
which indicate that the intrinsic size is resolved well be- 
yond the expected level of interstellar image broadening 
at A < 3.5 mm [3] and that the intrinsic size is smaller 
at higher frequencies. 

As mentioned in the Introduction, existing size mea- 
surements are subject to large uncertainties. While our 
best- fit profiles disfavor existing models, such models 
cannot be ruled out based on current data. Future, more 
accurate, measurements will allow us to draw more de- 
cisive conclusions. Our current analysis underlines the 
importance of future improvements in the size measure- 
ments, and provides a methodology for interpreting fu- 
ture results. 



IV. SUMMARY 

We have shown in § II. 2 that measurements of the spec- 
trum and of the wavelength-dependent size of the radio 
emission from Sgr A* indicate that this emission is dom- 
inated by optically-thick quasi-thermal plasma, and that 
the observed spectrum reflects the spatial dependence of 
the electron temperature (rather than the energy distri- 
bution of electrons at some particular radius). We have 
derived the electron temperature and magnetic field pro- 
files [Eqs. ([8|) and ([TTjl ]. and found that the electron tem- 
perature increases above the virial temperature beyond 
a distance of a few tens of Schwarzschild radii from the 
black hole. The observed rotation measure was then used 
to constrain the density profile [Eqs. ifTT)) and p8|) ] . The 
low density inferred for the gas near Sgr A* could in 
principle be accounted for by winds from the innermost 
S-stars [l6l | . Although we have not proposed a dynamical 
model for the plasma, we have pointed out in § II. 4 that 
observations are consistent with an isentropic gas profile 
and equipartition magnetic field. 

Our results imply that the radio emitting gas cannot 
be confined to a thin disk and the fiow geometry must be 
geometrically thick. Moreover, the radio-emitting gas is 
not likely to be infiowing but rather outflowing, since its 
thermal kinetic energy exceeds the gravitational binding 
energy beyond a radius of a few tens of Schwarzschild 
radii [ri3 > 2, see Table I and Eq. ([7])]. There may also 
be a colder accreting component that is sub-dominant 
in terms of its synchrotron emission. Such a component 
would likely be confined to a thinner disk geometry that 
would have only a limited effect on the rotation mea- 
sure of the radiation emitted by the hot outflowing at- 
mosphere above it. (However, if the cold component is 
optically-thick then it would make the image of SgrA* 
asymmetric at a level that would depend on the inclina- 
tion of the disk. The resulting frequency dependence of 
the image centroid location could be constrained by ob- 
servations.) An electron temperature profile which does 
not decline with increasing radius as fast as the virial 
temperature does, would be consistent with an adiabatic 
outfiow in which the electrons are hotter than the pro- 
tons because their relativistic temperature declines with 
decreasing density as T oc n^^^ with an adiabatic index 
(r = 4/3) that is smaller than that of the protons (5/3). 
Heat conduction could also help to flatten the electron 
temperature profile [20j . 

The large uncertainty in the inferred source size, re- 
flected by the large uncertainty in the index a = 1 ± 0.3 
of the relation r oc j/", translates to a large uncertainty 
in the temperature and magnetic field profiles, as implied 
by Eqs. pT|) and Our preliminary conclusions from 
existing data differ from current theoretical models [2^, 
although not at a statistically robust level. These po- 
tential discrepancies provide added incentive to obtain 
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better data through future observations. 

Finally, we note that the uncertainty in the determi- 
nation of the density profiles is related not only to un- 
certainties in the source size measurements, but also to 
the lack of simultaneous multi-frequency measurements 
of the rotation measure (§ II. 3). An accurate determi- 
nation of the density profile would require therefore not 
only accurate size measurements, but also simultaneous 
multi-frequency measurements of the rotation measure. 
Ultimately, direct imaging of Sgr A* with a Very Large 
Baseline Array at sub-mm wavelengths [^, [3] would 
resolve the accretion fiow near the black hole ISCO and 
unravel unambiguously the properties of the emitting gas 
there. 
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